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I. INTRODUCTION 

Quantum information processing devices are expected 
to be efficient tool for solving some practical problems 
which are exponentially hard for classical computers^. 
Their potential computational performance is achieved 
by exploiting quantum evolution of many particle system 
in exponentially large Hilbert space, necessarily including 
evolution steps through entangled states. Experimental 
implementation of Shor's quantum factoring algorithm in 
seven spin- 1/2 nuclei molecule have been demonstrated 2 . 

The question of whether a scalable implementation 
of quantum computer is possible in near future im- 
plies therefore the question of whether one can pro- 
tect the fragile entangled states from destructive envi- 
ronment. The dynamics of coherence loss of entangled 
many-particle clusters has attracted much attention re- 
cently. Some authors simulated the noisy environment 
as a single bosonic bath embracing whole dustes^^^. 
An alternative approach in which the noise sources act- 
ing on each cluster constituent are uncorrelated was also 
studied^. The realistic model of environment will be 
somewhere between these two cases. Still, the quanti- 
tative account for partially correlated environment com- 
plicates analysis much£, even for two particle system 9 . 
Until recently experimental data on decoherence of large 
clusters of highly entangled particles were also unavail- 
able. In 2004 the coherence dynamics of groups of up 
to 650 entangled nuclear spins was observed for the first 
time 10 . This paper is motivated by this experimental 
breakthrough indicating the partial correlation of the en- 
vironment. 

In this paper, we derive the dependence of decoherence 
rate of large spin clusters due to completely correlated 
and uncorrelated perturbation. The results are gener- 
alized to the system consisting of nuclear spins 1=1/2 
experimentally studied in the papesA^ by using solid-state 
NMR technique for powdered adamantane samples. Our 
results show that its dynamics resembles the decoherence 
due to interaction with a composite bath with a given ra- 
tio of correlated and uncorrelated terms. The dependence 
of decoherence rate on number of spins in the cluster was 



obtained. 

This paper is organized as follows. The investigated 
system and experimental procedures are described in Sec. 
ITT1 In Sec. IIIII we calculate the dynamics of NMR sig- 
nal for the cases of totally correlated/uncorrelated ex- 
ternal perturbations and for the experimental situation 
when decay is caused by internal dipole-dipole interac- 
tion. Comparison with experiment data and discussions 
are given m Sec. Concluding remarks are summa- 

rized in Sec. [V] 



II. SYSTEM 

Our study was stimulated by the recent results of Kro- 
janski and Sutesifi. In their experiments the system of 
nuclear spins-1/2 (protons) of the powdered admaman- 
tane sample was explored by methods of NMR. Initially 
a system, placed in the external magnetic field Hq along 
z axes, is in thermal equilibrium 
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where N is number of spins, 7 is spin gyromagnetic ratio, 
k is Boltzmann constant, T is temperature and P z is z 
component of i-th spin operator. With the help of spe- 
cial sequence of radio-frequency pulsesiS the high-order 




FIG. 1: (Color online). Sketch of random array of nuclear 
spins. Two clusters in partially correlated environments are 
shown. 
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correlations between spins grow thereby creating an en- 
semble of weakly coupled spin clusters. As a result, to 
describe evolution of spins in the sample it suffices to 
consider only the dynamics of one such cluster with well 
defined number of spins ni&lliiA. 

Existence of high-order coherences in n-spin system 
can be formally described by presence of the off-diagonal 
elements pij of the spin density operator in any rep- 
resentation whose basis states can be characterized by 
the total quantum magnetic numbers: Mi\i) = I z \i), 
= Ix\3)- Following the notation used in multi- 
ple quantum NMR experiments^ we say that every off- 
diagonal density matrix element p^ represents the coher- 
ence of the order M where M — Mi—Mj. The number of 
coherences (different off-diagonal elements) of the order 
M in n-spin system at large n is given by 
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It is conventional to assume that after long pulse se- 
quence spins are prepared in the state described by the 
density operator p(0) with all even coherences excited 
with equal probabilityi2*ii. 

After the system is prepared in this high-correlated 
state it decays under dipole-dipole interaction given by 
the Hamiltonian 



H. 



dd 



(Wit 



(3) 



where d jk = ±fr 2 7 2 (I - 3cos 2 [6> ife ])/r| fc and r jk , 6 jk are 
corresponding absolute value and the angle with z direc- 
tion of the vector connecting j-th and fc-th spins. 
The system, evolving according to 



p(t) = cxp (-~H d dtj p(0) cxp ( jHddt 



(4) 



does not produce experimentally observable signal. To 
analyze the effect of dipole-dipole interaction, it un- 
dergoes conversion step by another sequence of radio- 
frequency pulses described in RefJi. During this 
step multiple-quantum coherences are converted back to 
single-quantum longitudinal magnetization. After apply- 
ing a resonant frequency 7r/2 pulse which converts the 
longitudinal magnetization into transverse one the re- 
sulting longitudinal magnetization can be determined by 
measuring the free induction decay. The free induction 
decay amplitude right after 7r/2 pulse is proportional to 



S(t) ex Tr [p(t)p(0)} , 



(5) 



where t is the time the system freely evolved under 
dipole-dipole Hamiltonian between the end of the prepa- 
ration step and the beginning of the conversion step. The 
experiment has to be repeated for sequence of decay times 
t to obtain the decay of coherence. The overall signal can 
be presented as a sum of contributions corresponding to 



different coherence orders M— 



(6) 
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The decay times for Sm if) were also measured experi- 
mentally^ as a function of coherence order M for differ- 
ent cluster sizes n. 



III. THEORY 

A. Decay of NMR signal due to 
uncorrelated/correlated external baths 

First of all, consider a model when the decay of co- 
herence occurs due to interaction with the external bath. 
We do not specify the bath itself and use the generic pic- 
ture. In other words, in this subsection we consider the 
system without dipole-dipole interaction between spins. 
Instead we introduce some interaction with external bath 
which causes the initial spin coherence to decay. In this 
paper we focus only on the dephasing part of this inter- 
action. A classical analog of such model can be a cluster 
of spins in fluctuating external magnetic field directed 
along z axiA 

Henceforward, we use Zeeman basis |o) = \ai...a n ), 
where a, = ±1 and ZJ|cij) = (a,i/2)\a,i). If we consider the 
interaction of a single spin with the bath, its evolution in 
Zeeman basis is given by 

P±i,±iW = P±i,±i(0); px,-x{t) = pt^^e-r^, (7) 

where we used the interaction representation and the ex- 
plicit form of the decay function T(t) is determined by 
the nature of specific spin-bath interaction. Our main 
question is how the rate of the collective decoherence of 
the correlated spin cluster, measured by the technique 
given in the previous section, differs from the one for the 
single spin dynamics Q). The answer certainly depends 
on degree of correlation of the bath at different spin sites. 

As the first example, we consider the limiting case of 
completely uncorrelated environment: each spin interacts 
with its own bath assuming no correlations between baths 
related to different spins. In this case the matrix elements 
of n spin system density operator evolve according to 
Ref£ 



Pab{t) = Pab{0) exp (-T ab (t)) , 



(8) 



where collective decay function T a b(t) can be expressed in 
term of single spin decay function T(t) as r a t,(t) = fT(i) 
and / = (1/2) \ a i — &i| is t ne Hamming distance be- 
tween the spin states \a) and \b). The value of Hamming 
distance / has the same parity as coherence order M and 
is within the limits / G [M, n] . The number of configu- 
rations for given / and M for the system of n spin-1/2 
can be found as 
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We can calculate the observable decay of NMR signal 
S(t) according to I0I8I) as 

S(t) = ^2\p ab (Q)\ 2 exp[-r ab (t)}, (10) 

a . b 

where we need to carry out the summation over all pos- 
sible amplitudes \p a b\- The signal contributions Su(t) 
due to certain coherence order M can be evaluated by 
use of the same formula ((10(1 . Although in this case one 
needs to take the sum over only the subset of configu- 
rations {|a)(&|} € 9Jt for whose the additional condition 
^2j( a j — 2M is satisfied. The situation is greatly 

simplified by assuming that all even coherences are ini- 
tially excited with equal probability: |p a fc(0)| = const if 
(l/2)X}j( a j ~~ bj) = 0,2,4, ..; while all other coherences 
are not existent p a b(0) = 0i2iAi. We can write 

S M (t)<x ]T exp(-i/T(t)). (11) 

a,bCVJl 

Integrating over all / with corresponding weight @ we 
obtain, for nT(t) < 1, n > 1 and M < n/2, 

S M (t)=exp(-|r(t)). (12) 

We are interested in times up to 1/e decay time where 
formula (|12|l is valid. Moreover, since nT(t) < 1 the 
decay function is only in the onset regime: T(t) < 1/n <C 
1 for n ^> 1. Therefore, we take only the lowest non- 
vanishing order of decay function in time (which is always 
quadratic) 

r(t) = at 2 + 0(t 4 ). (13) 
Using 1)12113(1 we obtain 

S M (t) = cxp(-|crf 2 ) . (14) 

We emphasize that while we used the short-time expan- 
sion for the decay function 1(13(1 we expect the formula 
(|14|l to be valid up to 1/e decay time for the large num- 
ber of spins n> 1. 

As a second example, we consider the case of com- 
pletely correlated environment when the whole cluster 
interacts with the same external bath. The dynamics of 
the density matrix elements is given by£ 

Pab(t) = Pab(0) exp (~M 2 T(t)) , (15) 

and the signal Sjw(^) decays as 

S M (t) = cxp (-M 2 r(i)) ~ exp (-M 2 at 2 ) . (16) 

Formula 114(1 should be compared with 1(16(1 . Both re- 
sults show that the decay of signal SW (*) can be approx- 
imated by the Gaussian function up to 1/e decay times 
for n 3> 1. However, M dependencies for two formulas 
are totally different. The decay of Su(t) for uncorre- 
cted environment does not demonstrate any dependence 



on coherence order M while for correlated case it strongly 
depends on it. Thus, we established the distinctive fea- 
tures of the influence of correlated/uncorrelated environ- 
ments onto spin cluster dynamics which can be observed 
experimentally by NMR methods. 

B. Decay of NMR signal due to internal 
dipole-dipole interaction 

In the experiment by Krojanski and Suter— the deco- 
herence is caused not by external bath but due to integral 
dipole-dipole interaction between spins. However, as we 
show below the resulting behavior of the system can be 
interpreted with the help of results obtained in the pre- 
vious subsection. 

The dipole-dipole Hamiltonian 10 commutes with Zee- 
man Hamiltonian 

H z = - 7 hH J2 I i- (17) 
3 

However, the complexity of the system and especially 
the fact that two terms P Z I Z and F ■ I k do not com- 
mute makes impossible to find the exact (analytical or 
numerical) solution of the problemiii&ii. Existence of 
high order coherences in the state described by prepared 
density operator p(0) also complicates the application of 
the traditional method of moments which enables to de- 
scribe the decay of coherence without solving explicitly 
for eigenvalues and eigenstates of energy in case of single- 
quantum NMR experiments 14,15 . In our case the decay 
of signal is not proportional to the autocorrelation func- 
tion Tr{I x (t)I x }, as in the case of decay of free induc- 
tion signal^, but is given by density operator correla- 
tor (jSJ) where one needs to evaluate the summation of 
exponentially large number of terms. In order to ob- 
tain the analytical results, we focus on pure dephasing 
effect of dipole-dipole interaction neglecting any spin ex- 
change between spins that is described by flip-flop term 
I 3 X I X + I 3 y Iy in Hamiltonian Q. The dipole-dipole de- 
phasing Hamiltonian has the form 

H* dd = 2j2d jk P z I k z . (18) 

3<k 

Because dephasing is not associated with energy transfer 
mechanism it is generally the fastest source of decoher- 
enc o 16 i 17 . It becomes the sole process for decoherence 
in the limit of " unlike spins"^*^ when spin exchange is 
suppressed. The consideration of only this type of inter- 
action enables analytical calculations which are also jus- 
tified by good agreement with experiment in wide range 
of parameters as it will be demonstrated below. 

In the Zeeman representation the off-diagonal den- 
sity matrix elements evolve according to (JSJ where decay 
function is given by 

F a b(t) = ^X!^ fe ( a J afe - b J b k)- (19) 

j<k 
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momenfei* and degree of correlation p is denned as 
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FIG. 2: (Color online). Decoherence rate as function of coher- 
ence order for different spin cluster sizes. The points represent 
experimental values-^. The solid lines are obtained with ac- 
cordance to theoretical formula iYI'-it . Degree of correlation 
p and Van Vleck second moment evaluated from comparison 
with experiment are given in Table 1. 



The dynamics of normalized NMR signal Hl()(l can be 
analytically expressed as 



S(t) = 16]J J2 |Pa b (0)| 2 cos 2 (id fe ,i). 



(20) 



j<k a l ,b l =±l; 



Exact analytical expression l|2U|) does not provides us 
with much information yet. Specifically, we intend to ob- 
tain explicit dependence on number of spins in the clus- 
ter. For this purpose we again assume that all even co- 
herences arc initially excited with equal probability and 
the size of the cluster is large n 3> liSiii. After perform- 
ing some algebra (the details are given in Appendix) we 
obtain the expression for normalized signal 



n at 



S M (t) = 1- P M 2 — - (1 -P)~ + 0(t% (21) 

in the second order in time. Here a = M2/9 where M2 = 
(9/4:)h~ 2 ^ d 2 k is Van Vleck expression for the second 



TABLE I: Degree of correlation and second moment for 
CioHis obtained from decoherence rates for different cluster 
sizes. 



n 


26 


41 


71 


116 


189 


309 


477 


650 


M 2 , lOV 2 


1.50 


1.65 


1.60 


1.60 


1.65 


1.60 


1.60 


1.55 
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0.27 


0.28 


0.33 


0.33 


0.32 


0.32 


0.32 


0.32 
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d 2 

a jk 



(22) 



so that < p < 1. Formula (|21|l is valid only at short 
time scales nat 2 <C 1, while we are also interested in 
much larger times up to nat 2 ~ 1. However, expansion of 
the signal 5m(£) in higher orders in time becomes exceed- 
ingly difficult. Therefore, to continue l|21l) to the longer 
times we use the analogy with the investigated limiting 
cases (|14l 11611 . Formula (|21|) contains two terms propor- 
tional to M 2 and n/2 which can be regarded as contri- 
butions from correlated and uncorrelated perturbations 
to spin dynamics, respectively. In fact, the interaction 
described by Hamiltonian 118|) can be semiclassically in- 
terpreted as the perturbing magnetic field at the site of 
each spin (parallel or antiparallel to the strong external 
magnetic field) produced by all other spins in a cluster. 
The consequent spread of Larmor frequencies for differ- 
ent spins in the cluster causes destructive interference, or 
dephasing, observable by the decay of NMR signal. The 
limit of totally correlated perturbation p — 1 corresponds 
to the case djk = const leading to the same perturbing 
field for each spin in the cluster. In contrast, the case 
of absolutely random coefficients (djk)j = gives p = 
and fully uncorrelated dynamics. The realistic situation 
is expected to be in between these two limiting cases. 
Thus, we write (|2*T)l as 



S M (t) = pexp (-M 2 erf 2 ) + (1 -p)exp 



-at' 



(23) 



which is mathematically exact in up to the second order 
in time but continued to the longer times nat 2 ~ 1 . The 
total magnetic resonance signal from the cluster S(t) can 
be obtained by summation over all contributions from 
different coherence orders 5m (i) according to formulas 
© and lO 



S(t) 



P 



\Jnat 2 



= + (l-p)exp [-^af 



(24) 



In order to understand whether the obtained formulas 
()23l I24|l adequately describe the real experimental situ- 
ation we should check them with experiment data. The 
comparison of presented theory and experiment is given 
in the next section. 



IV. COMPARISON OF THEORY WITH 
EXPERIMENT AND DISCUSSIONS 

Recent experiment*^ allowed us to estimate the de- 
gree of correlation parameter for spin clusters in adaman- 
tane samples. In Fig. |3 we show curves of decay rates 
of various coherence orders for different cluster sizes fit- 
ted to experimental points. The decoherence rate was 
defined as the inverse of 1/e decay time and was eval- 
uated by solving the algebraic equation 5m(*) = l/ e 
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FIG. 3: (Color online). Decay of coherence from high- 
correlated spin clusters for different spin cluster sizes. The 
points represent experimental values—. The solid lines are 
values predicted by formula (12 It . 



where Sm (t) is given by formula The degree of 

correlation p and Van Vleck second moment Mi were ex- 
tracted with the use of MATLAB software by weighted 
least squares fitting to experimental data for every clus- 
ter size n. We minimized — /«) 2 /^? where Xi 
and /, are experimental points, f(xi) are corresponding 
theoretical solutions and Aj are experimental errors de- 
noted by vertical bars in Fig. [21 Obtained values of p 
and Mi are given in Table 1. As it follows from formula 
the definition of second moment it is determined by 
geometrical configurations and do not depend on cluster 
size n. Its moderate fluctuations around average value 
(If 2 = (1.60±0.05) • 10 9 s~ 2 ) can be attributed to experi- 
mental errors and corrections at small n. Obtained values 
for the second moment are comparable but not identical 
with the previous theoretical estimates and experimental 
measured values Mi ~ 2.6 • 10 9 s~ 2 for powdered solid 
adamant an a 1 8 1 1 9 . The difference can be explained either 
by crudity of the chosen model and neglecting flip-flop 
terms or by discrepancy in the adamantane samples used 
in different experiments. The question could be resolved 
by additional measurement of the second moment Mi for 
the given sample. 

Taking the average value of Mi = 1.6 • 10 9 s -2 and val- 
ues for p from Table 1 it is possible to predict the tempo- 
ral dependence of total NMR signal from high-correlated 
spin cluster l|24[l which was measured independently^ for 
different cluster sizes. The results shown in Fig. [3| are in 
good agreement with experiment. As can be seen from 
Fig. |3 the formula i|24|) describes the initial fast drop 



FIG. 4: (Color online). Examples of temporal dependence of 
the signal from high correlated spin cluster with size n=116 
and three values of degree of correlation p for perturbation: 
p=0 (dashed line, uncorrelated perturbation), 0.33 (solid line, 
partial correlation corresponding to experimental situation), 
1 (dotted line, correlated perturbation). The inset shows the 
decoherence rate as a function of coherence order M. 



of coherence with reasonable accuracy. The divergence at 
large times between formula (12411 (exact up second order 
in time) and experimental results can be attributed by 
the contribution of higher order terms. 

Let us note, that the measured values of decoherence 
rates for Sjw(£) give us only one time point in the tem- 
poral dynamics of the signal for each value of M and n. 
The parameters Mi and p obtained from fitting the so- 
lution of Sm (t) — 1/e to the experimental data allow us 
to reconstruct the total decay dynamics of Sm (t) up to 
1/e decay times. The following integration over all M 
provides the decay of overall signal S(t) which was mea- 
sured independently. This procedure does not, by any 
means, automatically guarantee the agreement of calcu- 
lated values of S(t) with experimentally measures ones. 
The correspondence of theoretical to experimental values 
demonstrates the good degree of consistency of presented 
theory. 

Formula Q24|l allows us to analyze the influence of de- 
gree of correlation on spin dynamics. Fig. 0] shows the 
decay of NMR signal for the spin cluster size of inter- 
mediate size n = 116 and three representative examples 
of degree of correlation p: p — (uncorrelated dynam- 
ics), p — 0.33 (partially correlated dynamics correspond- 
ing to the experimental situation) and p — 1 (correlated 
dynamics). One can see that initially all three curves 
decay equally. However, at later times the signal from 
the spin cluster subject to correlated perturbation ex- 
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hibits slower decay compared to uncorrelated perturba- 
tion. That result comes from the behavior of decoherence 
rate as function of coherence order M. As can be seen 
from inset of Fig.0] for uncorrelated perturbation all co- 
herence orders decay with the same, comparatively high, 
rate (na/2) 1 / 2 . In contrast, the decay rate for correlated 
spin dynamics increases linearly with absolute value of 
M as a 1 / 2 !-/^. For the most probable configurations, 
which according to © are those with M w 0, the decay 
rate for correlated perturbation is actually less than that 
for uncorrelated perturbation. The fact that correlated 
environment is acting more delicate on specific groups 
of states is not surprising. In particular, quantum com- 
puting error avoiding schemes based on decoherence free 
subspaces^SH are based on this property. 

For implementation of large-scale quantum computa- 
tion the scaling of decoherence rate with number of qubits 
is important. From the expression 1)24(1 it transpires that 
decoherence rate of a spin cluster defined as inverse 1/e 
decay time always increases as oc \fn with number of 
spins n although the corresponding factor depends on 
degree of correlation p. The square root of n scaling was 
indeed experimentally discovered recently by Krojanski 
and Suter—. 

For quantum information processing applications it is 
also important to evaluate the error of a quantum com- 
puter, represented by a cluster of high correlated spins, 
induced by dipole-dipole interaction between spins. The 
error is defined as deviation of NMR signal from its ini- 
tial value due to decoherence processes during the time 
required for elementary gate operation t g : S n = 1 — S(t g ). 
In order to provide successful implementation of quantum 
error correction schemes, one needs to maintain this er- 
ror below the small threshold guarantying fault-tolerance 
operation of these procedures^. Taking the smallness of 
the parameter S n into account one can use l(24() to obtain 

S n oc nt 2 . (25) 

This shows that if the error is small it scales linearly with 
number of spins independently of degree of correlation. 
The linear scaling of error agrees with theoretical results 
for bosonic models of environment^ and suggests that 
the worst case scenario of " superdecoherence"— is not re- 
alized for this particular system. 

V. SUMMARY 

In summary, we have presented a theory of coherence 
decay of entangled spin clusters states due to internal 



dipole-dipole interactions. Its dynamics resembles the 
decoherence due to interaction with a composite bath 
consisting of fully correlated and uncorrelated parts. The 
perturbation due to correlated terms leads to the slower 
decay of coherence at larger times. The decoherence rate 
scales up as a square root of the number of spins giving 
the linear scaling of the resulting error. The results ob- 
tained can be useful in analysis of decoherence effects in 
spin-based quantum computers. 
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APPENDIX A: SDF 

The signal contributions SmC^) can be evaluated by 
use of the formula 1(1 Ojl 

S M (t)= IPafc(0)| 2 ex P [-r Qb (i)] (Al) 

Here domain 971 denote all configurations \a)(b\ = 
|oi . .. a n )(bi . . .b n \ related to the certain coherence or- 
der M: 

J2(a j ~b j )=2M. (A2) 

3 

For every configuration \a)(b\ we can divide the total set 
of n spins into two subsets E and N, 

en = h, Vi e E; aj = -b h Vj G N. (A3) 

By the use of definitions l(A3(l the decay function T a i,(i) 
can be simplified as 
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T ab (t) = i- djkjajak - bjb k ) = i- ^ d 3 k{a 3 a k - b 3 b k ) + ij ^ d jk (a :j a k - bjb k ) 



j,k 



+ i 



I E d i k> 



jeE,keE 

a 3 a k - bjbk) = H djkdjak. 



j£N,k£f 



(A4) 



r 



Assuming that all even coherences are initially excited 
with equal probability, namely |/?ab(0)| = const for all 
s,i e 3)1 (and for all other even order coherences), we 
obtain the following formula for the decay of Sm (t) ac- 
cording to HA1IA4|) 

S M (t) oc 2J cxp[-i< 2J d 3k a 3 a k ]. (A5) 

We can redistribute the summation in (|A5|I in the follow- 
ing way 



E =E E • 

a,6cOT E.N a 3 ,a fc 



(A6) 



Here the first sum in the left part of the equation is over 
all possible choices of subsets E, N in the set of n spins 
and second sum is over all possible values of aj , a k for 
i e E, k e N. It is easy to see that values aj for j G E can 
take any values aj — ±1 since they do not contribute to 
(IA2|) and, therefore, do not change coherence order M. 
We still have the condition for the values a k , k £ N: 



(A7) 



kef 



Thus, we can evaluate the summation over aj for jeE 
first and obtain 

s M (t) « E E Il cos (*E d ^ afc ) ( A8 ) 

E,N a k ,kefijeE keti 

« E E (i-yE(E^) 2 + °^ 4 )V 

E,Na fc ,fc£N y jeE fcGN J 

Now we consider the term ( S dj k a k ) 2 to express it 

j'eE fceN 

in terms of parameters of the material. We write 



dj k = dj + Sjk, 



(A9) 



where average coupling constant is defined as 



jk- 



(A10) 



ken 



Here / is Hamming distance between \a) and \b) or the 
number of spins in subset N. Note, that X^feeN djk = 0- 
By use of (|A7IA9|I we obtain 

d^) 2 = (djfM 2 + (All) 
fcGN ken 
where we neglected cross-terms X^fceN dj k a k whose contri- 
bution is negligible for the large cluster sizes. For n ^S> 1 
we also approximate the summation over subsets by the 
summation over total cluster with correction to the num- 
ber of terms in the sum: £N — (n/(n — /)) X^"eE = 
( n / f)J2jen- These assumptions should lead to the 
asymptotically correct value of Suit) for n — ► oo. We 
then evaluate the the signal decay as 



S M {t) oc 2J 



e.:<i 



1 - 



t 2 n-f 



( P M 2 + (1-P)/)E^ 



-o(t 4 ), 



where the parameter p is defined as 

^4(E>) 2 /E/: 



(A12) 



(A13) 



After integration over all possible / we deduce the closed, 
analytical form for the signal, exact up to second order 
in time. 



Su(t) = 1 - - ^ (pA/ 2 + (1 - p)|) + 0(t 4 ), (A14) 

Here M2 — (9/4)fi, -2 JN d 2 k is Van Vleck expression for 
the second momen1ji£. 
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